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Abstract 

Motivated by forthcoming experiments at RHIC and LHC, and results from SPS, a 
review is given of the present state of event-by-event fluctuations in ultrarelativistic 
heavy-ion collisions. Fluctuations in particle multiplicities, ratios, transverse mo- 
menta, rapidity, etc. are calculated in participant nucleon as well as thermal models. 
The physical observables, including multiplicity, kaon to pion ratios, and transverse 
momenta agree well with recent NA49 data at the SPS, and indicate that such stud- 
ies do not yet reveal the presence of new physics. Predictions for RHIC and LHC 
energies are given. The centrality dependence with and without a phase transition 
to a quark-gluon plasma is discussed - in particular, how the physical quantities 
are expected to display a qualitative different behavior in case of a phase transi- 
tion, and can be signaled by anomalous fluctuations and correlations in a number 
of observables. 
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1 Introduction 

The importance of event-by-event physics is evident from the following simple 
analogy: Stick a sheet of paper out of your window on a rainy day. Keeping 
it there for a long time - corresponding to averaging - the paper will become 
uniformly wet and one would conclude that rain is a uniform mist. If, however, 
one keeps the sheet of paper in the rain for a few seconds only, one observes 
the striking droplet structure of rainQ. Incidentally, one has also demonstrated 
the liquid-gas phase transition! Analyzing many events gives good statistics 
and may reveal rare events as snow or hail and thus other phase transitions. 
The statistics of droplet sizes will also tell something about the fragmentation, 
surface tension, etc. By varying initial conditions as timing and orienting the 
paper, one can further determine the speed and direction of the rain drops. 

Central ultrarelativistic collisions at RHIC and LHC are expected to produce 
about ~ 10 4 particles, and thus present one with the remarkable opportunity 
to analyze, on an event-by-event basis, fluctuations in physical observables 
such particle multiplicities, transverse momenta, correlations and ratios. Anal- 
ysis of single events with large statistics can reveal very different physics than 
studying averages over a large statistical sample of events. The use of Hanbury 
Brown-Twiss correlations to extract the system geometry is a familiar appli- 
cation of event-by-event physics in nuclear collisions [1] and elsewhere, e.g, in 
sonoluminoscence [2]. The power of this tool has been strikingly illustrated 
in study of interference between Bose-Einstein condensates in trapped atomic 
systems [3]. Fluctuations in the microwave background radiation as recently 
measured by COBE [4] restrict cosmological parameters for the single Big 
Bang event of our Universe. Large neutrons stars velocities have been mea- 
sured recently [5] which indicate that the supernova collapse is very asymmet- 
rical and leads to large event-by-event fluctuations in "kick" velocities during 
formation of neutron stars. 

The tools applied to study these phenomena do, however, vary in order to op- 
timize the analysis and due to limited statistics. The COBE and the interfer- 
ence in Bose-Einstein condensates require study of fluctuations within a single 
event. The HBT studies in heavy- ion collisions and sonoluminosence requires 
further averaging over many events in order to obtain sufficient statistics; one 
has not yet studied fluctuations in source radii event-by-event. Anisotropic flow 
requires an event-plane reconstruction in each event [6] but again averaging 
over many events is necessary to obtain a statistically relevant measurement 
of the flow. The event-by-event fluctuations in heavy-ion collisions (and neu- 
tron star kick velocities) go a step further by studying variations from event 
to event. 



Originally, this analogy was given by Prof. A.D. Jackson 
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Studying event-by-event fluctuations in ultrarelativistic heavy ion collisions 
to extract new physics was proposed in a series of papers [7-9], in which the 
analysis of transverse energy fluctuations in central collisions [10] was used to 
extract evidence within the binary collision picture for color, or cross-section, 
fluctuations. More recent theoretical papers have focussed on different aspects 
of these fluctuations, such as searching for evidence for thermalization [11-13], 
correlations between transverse momentum and multiplicity [14], critical fluc- 
tuations at the QCD phase transition [15-19] and other correlations between 
collective quantities [20]. 

Intermittency [21] studies of factorial moments of multiplicities are related to 
event-by-event fluctuations. One of the motivations for intermittency studies 
was the idea of self-similarity on small scales, an idea borrowed from chaos 
theories. The factorial moments of particle multiplicities did find approximate 
power law behavior when the intervals of rapidity and angles were made in- 
creasingly smaller, at least until a certain small scale. The power law scaling 
in nucleus-nucleus collisions was, however, weaker than in proton-proton colli- 
sions. This indicated that the stronger correlations in proton-proton collisions 
were mainly due to resonances, minijets and other short range correlations, 
but that they were averaged out in nuclear collisions by summing over the 
many individual participating nucleons. The scaling was not a collective phe- 
nomenon and indications of new physics were not found [22]. In more recent 
event-by-event fluctuation studies the self-similar scaling idea is abandoned. 
They concentrate on the mean and the variance of the particle multiplicities 
per event and correlations between different particle species, transverse mo- 
mentum, azimuthal angle, etc. One directly compares to expectations from 
proton-proton collisions scaled up by the number of participants. One fol- 
lows these fluctuations and correlations for heavy-ion collisions as function of 
centrality and system size searching for anomalous behavior as compared to 
proton-proton collisions. 

Recently NA49 has presented a prototypical event-by-event analysis of fluctu- 
ations in central Pb+Pb collisions at 158 GeV per nucleon at the SPS, which 
produce more than a thousand particles per event [11]. The analysis has been 
carried out on ~100.000 such events measuring fluctuations in multiplicities, 
particle ratios, transverse momentum, etc. 

Results from the RHIC collider are eagerly awaited [23]. The hope is to observe 
the phase transition to quark-gluon plasma, the chirally restored hadronic mat- 
ter and/or deconfinement. This may be by distinct signals of enhanced rapid- 
ity and multiplicity fluctuations [24,17] in conjunction with J/\l/ suppression, 
strangeness enhancement, rj' enhancement, constant (critical) temperatures 
vs. transverse enery or rapidity density [25], transverse flow or other collective 
quantities as function of centrality, transverse energy or multiplicity as will be 
discussed in detail below. 
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The purpose of this review is to understand these and other possible fluctu- 
ations. We find that the physical observables, including multiplicity, kaon to 
pion ratios, and transverse momenta agree well with recent NA49 data at the 
SPS, and indicate that such studies do not yet reveal the presence of new 
physics. Predictions for RHIC and LHC energies are given. The centrality 
dependence with and without a phase transition to a quark-gluon plasma is 
discussed - in particular, how the physical quantities are expected to display 
a qualitative different behavior in case of a phase transition, and how a first 
order phase transition could be signaled by very large fluctuations. 



2 Phase Transitions and Fluctuations 

Lattice QCD calculations find a phase transition in strongly interacting matter 
which is accompanied by a strong increase of the number of effective degrees 
of freedom [26,27]. The Early Universe underwent this transition at a time 
t = 0.3 — 0A(T c /MeV)~ 2 seconds. For a hadronic gas melting temperature 
of T c = 150 MeV this occurred around 15 microseconds after the Big Bang. 
By colliding heavy nuclei we expect to reproduce this transition at sufficiently 
high collisions energies. 



2.1 Order of the QCD phase transition 

The nature and order of the transition is not known very well. Lattice calcu- 
lations can be performed for zero quark and baryon chemical potential only, 
Hb = 0, where they suggest that QCD has a weak first order transition pro- 
vided that the strange quark is sufficiently light [26,27], that is for 3 or more 
massless quark flavors. The transition is due to chiral symmetry restoration 
and occur at a critical temperature Tc — 150 MeV. In pure SU(3) gauge the- 
ory (that is no quarks, Nf = 0) the transition is a deconfinement transition 
which is of first order and occurs at a higher temperature T c ~ 260 MeV. 

However, when the strange or the up and down quark masses become massive, 
the QCD transition changes to a smooth cross over. The phase diagram is 
then like the liquid-gas phase diagram with a critical point above which the 
transition goes continuously through the vapor phase. For reasonable values 
for the strange quark mass, m s ~ 150 MeV and small up and down quark 
masses, lattice calculations find either a weak first order transition [26,27] 
or a smooth soft cross-over [28]. In case of a weak first order transition, the 
latent heat and density discontinuities and the signals, that depend on these 
quantities, will be small. 
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For exactly two massless flavors, m U) d = and m s = oo, the transition is 
second order at small baryon chemical potential. Random matrix theory finds 
a 2nd order phase transition at high temperatures which, however, change into 
a 1st order transition above a certain baryon density - the tricritical point. For 
small up and down quark masses the transition changes to a continuous cross- 
over at zero baryon chemical potential but remains a first order at large baryon 
chemical potential. A critical point must therefore exist at small but finite 
baryon chemical potential which may be searched for in relativistic heavy-ion 
collisions [17-19]. 



T 



Early Universe 

LHC: Pb-Pb (2005) 
(s = 5.5 ATeV 



Quark-Gluon Plasma 
deconfinement 



RHIC: Au-Au (now!) 
{s = 200 AGeV 

SPS: Pb-Pb (now) 
E b = 158 AGeV 
(s= 17 AGeV 



E b = 40 AGeV 




Hadron Gas 
confinement 



AGS: Au-Au, (past) 
E b = 11 AGeV 
fs = 5 AGeV 



neutron stars? 



Fig. 1. An illustration of the QCD phase diagram, temperature vs. baryon chemical 
potential. The regions of the phase diagram probed by various high energy nuclear 
collisions are sketched by arrows. From [29]. 



2.2 Density, rapidity, temperature and other fluctuations 



Fluctuations are very sensitive to the nature of the transition. In case of a 
second order phase transition the specific heat diverges, and this has been 
argued to reduce the fluctuations drastically if the matter freezes out at T c [15- 
18]. For example, the temperature fluctuations have a probability distribution 
[30] 

AT 

w ~ exp(-CV(— ) 2 ) ; (1) 
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Fig. 2. Sketch of droplet formation (top) in a continuous background of hadron vs. 
rj. The corresponding rapidity distributions (bottom) are shown for the continuous 
hadronic background with and without droplets (target and projectile fragmentation 
regions are excluded). 



a diverging specific heat near a 2nd order phase transition would then remove 
fluctuations if matter is in global thermal equilibrium. The implications of such 
critical phenomena near second order phase transitions and critical points are 
discussed in detail in [17,18]. It is found that the expansion of the systems slows 
the growth of the correlation lengths associated with the critical phenomena 
and the systems "slows out of equilibrium", which affects the experimental 
signatures related to transverse momenta and temperatures. 

First order phase transitions are contrarily expected to lead to large fluctua- 
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tions due to droplet formation [31] or more generally density or temperature 
fluctuations. These hot droplets will expand and hadronize in contrast to cold 
static quark matter droplets that may exist in cores of neutron stars [32]. 
In case of a first order phase transition relativistic heavy ion collisions lead 
to interesting scenarios in which matter is compressed, heated and undergoes 
chiral restoration. If the subsequent expansion is sufficiently rapid, matter will 
pass the phase coexistence curve with little effect and supercool [33,34]. This 
suggests the possible formation of "droplets" of supercooled chiral symmetric 
matter with relatively high baryon and energy densities in a background of low 
density broken symmetry matter. These droplets can persist until the system 
reachs the spinodal line and then return rapidly to the now-unique broken 
symmetry minimum. A large mismatch in density and energy density seems 
to be a robust prediction for a first order transition at large baryon densities. 
At high temperatures, which is more relevant for relativistic heavy-ion colli- 
sions (see Fig. (1)), the transition is probably at most weakly first order as 
discussed above. 

Density fluctuations may appear both for a first order phase transition and 
for a smooth cross-over. If the transitions is first order, matter may supercool 
and subsequently create fluctuations in a number of quantities. Density fluc- 
tuations in the form of hot spots or droplets of dense matter with hadronic 
gas in between is a likely outcome (see Fig. (2)). Even if the transition is a 
smooth cross-over, the resulting soft equation of state has a small sound speed, 
c 2 s = dP/de. The equation of state P(e) has in both cases a flat region that 
may be hard to distinguish in a finite systems existing for a short time only. We 
do not know the early non-equilibrium stages of relativistic nuclear collisions 
and the resulting initial density fluctuations, hot spots, etc. If the system be- 
comes thermalized at some stage, then a smaller c 2 s is likely to allow for larger 
density fluctuations since the pressure difference is smaller. Furthermore, in 
the subsequent expansion the density fluctuations are not equilibrated as fast 
when Cg is small because the pressure differences, that drive the differential 
expansion, are small. The dissipation of an initial density fluctuations can 
be estimated by a stability analysis [35]. Linearizing the hydrodynamic equa- 
tions in small fluctuations around the Bjorken scaling solution, an entropy 
fluctuation is typically damped by a factor (see Appendix A for details) 



Here, a typical formation time is r ~ 1 fm/c and freezeout time Tf ~ 8 fm/c 
as extracted from HBT studies [1]. The eigenvalues \± depend on the sound 
speed and the wave length of the rapidity fluctuations. As described in more 
detail in Appendix A, one of the eigenvalues are small and vanish for c s = 0. 
The resulting suppression of an initial density fluctuations during expansion is 
typically smaller than a factor 0.5. If density fluctuations are enhanced initially 




(2) 
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due to a softening of the equation of state due to smooth cross over, then they 
will largely be preserved later on. Yet, such fluctuations will be smaller than 
for a true first order transition forming supercooled droplets. 

Let us assume that hadrons emerge from a collection of density fluctuations 
or droplets with a Boltzmann distribution with temperature T and from a 
more or less continous background obeying approximate Bjorken scaling. The 
resulting particle distribution is 

oc Y fi e~ mt c °My-Vi)/T + background . (3) 
dyd 2 p t Y 

Here, y is the particle rapidity and p t its transverse momentum, /, is the 
number of particles hadronizing from each droplet i, and 

1, U + Zi 1 1 + Vi 

^ = -2 l0g — = 2 l0g — V t (4) 

is the rapidity of droplet i. The size, number and separation between droplets 
or density fluctuations will depend on the violent initial conditions. Between 
droplets a relatively continuous background of hadrons is expected in coexis- 
tence. In (3) the droplet is assumed not to expand internally neither longitu- 
dinally nor transversely. If it does expand, the emerging hadrons will have a 
wider distribution of rapidities which will be harder to distinguish from the 
background. 

When rrit/T ^> 1, we can approximate cosh(y — rji) ~ 1 + |(y — rji) 2 in Eq.(3). 
The Boltzmann factor determines the width of the droplet rapidity distribu- 



tion as ~ \jT/m t . The rapidity distribution will display fluctuations in rapid- 
ity event by event when the droplets are separated by rapidities larger than 



\r)i — r]j\£jyT/m t . If they are evenly distributed by smaller rapidity differences, 
the resulting rapidity distribution (3) will appear flat. 

The droplets are separated in rapidity by 1 77^ — r]j\ ~ Az/tq, where Az is the 
correlation length in the dense and hot mixed phase and tq is the invariant 
time after collision at which the droplets form. Assuming that Az ~ lfm — a 
typical hadronic scale — and that the droplets form very early To^lfm/c, we 
find that indeed — rjj\^\jT /m t even for the light pions. If strong transverse 
flow is present in the source, the droplets may also move in a transverse direc- 
tion. In that case the distribution in p t may be non-thermal and azimuthally 
asymmetric. 

Even if the transition is not first order, fluctuations may still occur in the 
matter that undergoes a transition. The fluctuations may be in density, chiral 
symmetry [36], strangeness, or other quantities and show up in the associated 
particle multiplicities. The "anomalous" fluctuations depend not only on the 
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Fig. 3. Charged particle density per participant at midrapidity are shown vs. cm. 
energy for pp and pp data [39-41] with open symbols. Full curve presents a linear fit 
to pp and pp data up to RHIC energies. Pb+Pb data from NA49 [11] and Au+Au 
data from PHOBOS [23] per participant are shown with filled symbols, and exceed 
the pp and pp by ca. 30% (dashed curve) . 

type and order of the transition, but also on the speed by which the collision 
zone goes through the transition, the degree of equilibrium, the subsequent 
hadronization process, the amount of rescatterings between hadronization and 
freezeout, etc. It may be that any sign of the transition is smeared out and 
erased before freezeout. That no anomalous event-by-event fluctuations have 
been found at CERN [24] within experimental accuracy indicate that no tran- 
sition took place or that the signals were erased before freezeout. Whether 
they remain at RHIC is yet to be discovered and we shall provide some tools 
for the analysis in the following sections. 



3 Multiplicity Fluctuations in Relativistic Heavy-Ion Collisions 



In order to be able to extract new physics associated with fluctuations, it is 
necessary to understand the role of expected statistical fluctuations. Our aim 
here is to study the sources of these fluctuations in collisions. As we shall see, 
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the current NA49 data (see Fig. (6)) can be essentially understood on the ba- 
sis of straightforward statistical arguments. Expected sources of fluctuations 
include impact parameter fluctuations, fluctuations in the number of primary 
collisions and in the results of such collisions, fluctuations in the relative ori- 
entation during the collision of deformed nuclei [10], effects of rescattering of 
secondaries, and QCD color fluctuations. Since fluctuations in collisions are 
sensitive to the amount of rescattering of secondaries taking place, we dis- 
cuss in detail two limiting cases, the participant or "wounded nucleon model" 
(WNM) [37], in which one assumes that particle production occurs in the in- 
dividual participant nucleons and rescattering of secondaries is ignored, and 
the thermal limit in which scatterings bring the system into local thermal 
equilibrium. 

Data at AGS, SPS and RHIC energies show that multiplicities are enhanced 
by ~30% in central collisions between heavy (A ~ 200) nuclei as compared 
to the WNM prediction (see Fig. (3)). Whether rescatterings increase relative 
fluctuations through greater production of multiplicity, transverse momenta, 
etc., or decrease fluctuations by involving a greater number of degrees of free- 
dom, is not immediately obvious [7,8]. Rescatterings probably increase both 
the average multiplicity and its variance but whether the relative amount 
of fluctuations are increased is model dependent. It has even been found in 
relativistic heavy ion collisions that the multiplicity fluctuations increase in 
the first few rescattering but then decrease again as the thermal limit is ap- 
proached. VENUS simulations [38] showed that rescattering had negligible 
effects on transverse energy fluctuations. 

We first review known multiplicities and fluctuations in the basic pp collisions, 
go on to study nucleus-nucleus collisions, and finally show in a simple model 
how phase transitions are capable of producing very significant fluctuations in 
particle multiplicities. 



3. 1 Charged particle production in pp and pp reactions 

Participant models or WNMs are basically a superposition of NN collisions. 
Such models have been studied extensively at these energies within the last 
decades at several particle accelerators and we here give a brief compilation 
of relevant results. 

The average number of charged particles produced in high energy pp and 
ultrarelativistic pp collisions can be parametrized by 




(5) 
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Fig. 4. The total number of charged particles produced in pp and pp collisions vs. 
cms energy s 1 / 2 . Data from bubble chamber [40], ISR [42], UA5 [41] and FNAL 
E735 [43]. 

for cms energies \fs^L,2 GeV. At ultrarelativistic energies the charged particle 
production is very similar in pp, pn and pp collisions and the parametrization 
of Eq. (5) applies in a wide range of cms energies 2 GeV^s 1 ^ 2 ^ 2 TeV as 
shown in Fig. (4). At SPS, RHIC and LHC energies, y/s ~ 20, 200, 5000 GeV, 
we find (Nch) — 7.3, 20, 60, respectively. 

At high energies KNO scaling [44] is a good approximation. KNO scaling im- 
plies that multiplicity distributions are invariant when scaled with the average 
multiplicity. Thus all moments scale like 



(6) 



at high energies where c q are constants independent of collision energy. The 
fluctuations, 



u N = «iV 2 > - (N) 2 )/(N) 



(7) 



therefore scale with average multiplicity, (N) , and therefore increase with col- 
lision energy as in Eq. (5). The fluctuations in the charged particle multiplicity 
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Fig. 5. The fluctuations in the total number of charged particles produced in pp and 
pp collisions from bubble chamber experiments [40], ISR [42], and UA5 [41]. Note 
the large difference between pp and thermal fluctuations uj 1 ^ ~ 2.2 at very high 
energies and the "accidental" crossing with pp fluctuations around SPS energies. 

can be parametrized rather accurately by 

u Nch ~ 0.35 — — (8) 

V*ch) 

as shown in Fig. (5) for pp and pp collisions in the same wide range of energies. 
At the very high energies breakdown of KNO scaling has been observed in the 
direction that the fluctuations are slightly larger. At SPS, RHIC and LHC 
energies we find uiN ch — 2.0, 6.2, 20, respectively in pp and pp collisions. 

In nuclear (AA) collisions the number of participating nucleons N p grow with 
centrality and nuclear mass number A. Therefore the average charged particle 
multiplicity and variance grows with N p , whereas the ratio and therefore the 
fluctuation to^ ch is independent of N p , and equal to the fluctuations in pp col- 
lisions. (Other fluctuations such as impact parameter will be included below.) 
Higher moments of the multiplicity distributions are large in high energy pp 
and pp collisions due to KNO scaling but in nuclear collisions such higher 
moments are suppressed by factors of 1/N p and are therefore less interesting 
than the second moment. This justifies our detailed analyzes of the variance 
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(or rms width) of the fluctuations. 



3.2 Fluctuations in the participant model 

In the participant or wounded nucleon models nucleus-nucleus collisions at 
high energies are just a superposition of nucleon- nucleon (NN) interactions. 
In peripheral collisions there are only few NN collisions, the collision zone is 
small, rescatterings few and the WNM should therefore apply. For central nu- 
clear collisions, however, multiple NN scatterings, energy degradation, rescat- 
terings between produced particles and other effects complicate the particle 
production and do enhance the multiplicities by ca. 30% as seen in experiment 
(see Fig. (3)). Thermal models may better describe central collisions as will be 
investigated afterwards. Yet, the WNM provides a simple baseline to compare 
to, when going from peripheral towards central collisions. 

Let us first calculate fluctuations in the participant model. Although the mul- 
tiplicities are somewhat underestimated, the measured multiplicity and trans- 
verse energy in nuclear collisions at AGS and SPS energies are known to scale 
approximately linearly with the number of participants [45,11]. In this picture 



where N p is the number of participants and is the number of particles 
produced in the acceptance by participant i. In the absence of correlations 
between N p and n, the average multiplicity is (N) = (N p )(n). For example, 
NA49 measures charged particles in the rapidity region 4 < y < 5.5 and finds 
(N) ~ 270 for central Pb+Pb collisions. Finite impact parameters (6^,3.5 fm) 
as well as surface diffuseness reduce the number of participants from the total 
number of nucleons 2A to (N p ) ~ 350 estimated from Glauber theory; thus 
(n) ~ 0.77. Squaring Eq. (9) and again assuming no correlations between 
different wounded nucleon emission, (jiiUj) = (ni)(rij) for % ^ j, we find the 
multiplicity fluctuations (see Appendix B) 

u N = u n + (n)u Np , (10) 




where oon, oo n and un are the multiplicity fluctuations in the total number of 
particles (within the acceptance), in each source, and in the number of sources 
respectively. 

A major source of multiplicity fluctuations per participant, u n , is the lim- 
ited acceptance. While each participant produces v charged particles, only a 
smaller fraction / = (n)/{u) are accepted. Without carrying out a detailed 
analysis of the acceptance, one can make a simple statistical estimate assum- 



(9) 
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ing that the particles are accepted randomly, in which case n is binomially 
distributed with cr(n) = vj(\ — /) for fixed v. Including fluctuations in v we 
obtain, similarly to Eq. (10), 



u n = 1 - / + fa* ■ 



(11) 



In nucleon-nucleon collisions at SPS energies, the charged particle multiplicity 
is ~ 7.3 and oo u ~ 1.9 [46]; as the multiplicity should be divided between the 
two colliding nucleons, we obtain (v) ~ 3.7 and thus / = (n)/(u) = 0.21 for 
the NA49 acceptance. Consequently, we find from Eq. (11) that uo n ~ 1.2. 
The random acceptance assumption can be improved by correcting for known 
rapidity correlations in charged particle production in pp collisions [42,40]. 

Multiplicities generally increase with centrality of the collision. We will use 
the term centrality as impact parameter b in the collision. It is not a directly 
measurable quantity but is closely correlated to the transverse energy pro- 
duced Et, the measured energy in the zero degree calorimeter and the total 
particle multiplicity N measured in some large rapidity interval. The latter is 
within the WNM approximately proportional to the number of participating 
nucleons 



For sharp sphere nuclei the number of participants drops from N p (b = 0) = 2 A 
in central collisions to N p (b = 2R) = in grazing collisions. For realistic nuclei 
with diffuse surface and with collision probabilities given by Glauber theory, 
the number of participants are 5-10% smaller in central collisions but slightly 
larger in peripheral collisions. 

As a consequence of nuclear correlations, which strongly reduce density fluc- 
tuations in the colliding nuclei, the fluctuations u>N p (b) in N p are very small 
for fixed impact parameter b [9]. Almost all nucleons in the nuclear overlap 
volume collide and participate. [By contrast, the fluctuations in the number 
of binary collisions are non-negligible.] Cross section fluctuations play a small 
role in the WNM [9]. Fluctuations in the number of participants can arise 
when the target nucleus is deformed, since the orientations of the deformation 
axes vary from event to event [47]. The fluctuations, u>n p , in the number of 
participants are dominated by the varying impact parameters selected by the 
experiment. In the NA49 experiment, for example, the zero degree calorime- 
ter selects the 5% most central collisions, corresponding to impact parameters 
smaller than a centrality cut on impact parameter, b c ~ 3.5 fm. We have 




overlap 



(12) 




o 



(13) 
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where (N p ) = (1/irb 2 ) J bc d 2 bN p (b) . The number of participants for a given 
centrality, calculated in [48], can be approximated by N p (b) ~ N P (0)(1 — b/2R) 
for < 6^,3.5 fm; thus 




For NA49 Pb+Pb collisions with N p (0) ~ 400 and (b c /2R) 2 ~ 5% we find 
u;jv p ~ 1.1. Impact parameter fluctuations are thus important even for the 
centrality trigger of NA49. Varying the centrality cut or b c to control such im- 
pact parameter fluctuations (14) should enable one to extract better any more 
interesting intrinsic fluctuations. Recent WA98 analyzes confirm that fluctu- 
ations in photons and pions grow approximately linearly with the centrality 
cut (b c /2R) 2 [49] as predicted by Eq. (14). The impact parameter fluctuations 
associated with the range of the centrality cut, such at total transverse en- 
ergy or multiplicity, can therefore be removed. However, fluctuations in impact 
parameter may still remain for a given centrality. The Gaussian multiplicity 
distribution found in central collisions changes for minimum bias to a plateau- 
like distribution [10]. 

Calculating ujn for the NA49 parameters, we find from Eq. (10), ujn ~ 1.2 + 
(0.77) (1. 1) = 2.0, in good agreement with experiment, which measures a mul- 
tiplicity distribution oc exp[-(iV - (N)) 2 /2(N)lu^ p ], where is 2.01 [11] 
(see Fig. (6)). 



3.3 Fluctuations in the thermal model 



Let us now consider, in the opposite limit of strong rescattering, fluctuations 
in thermal models. In a gas in equilibrium, the mean number of particles per 
bosonic mode n a is given by 

(n a ) = (exp(E a /T)-l)- 1 , (15) 

with fluctuations 

u na = l + (n a ) . (16) 
The total fluctuation in the multiplicity, N = J2a n a, is 

< = i + EWVEK). (i7) 

a a 

If the modes are taken to be momentum states, bosons/fermions have ther- 
mal fluctuations, ojn = 1 ± (n 2 ) / (n p ) where n p = (exp(e p /T) =)= 1) _1 is the 
boson/fermion distribution function, which are slightly larger /smaller than 
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Fig. 6. Event-by-event fluctuations of multiplicity (top) and pt (bottom) measured 
by NA49 in central Pb+Pb collisions at the SPS [11]. 
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those of Poisson statistics for a Boltzmann distribution, ojn = 1. The result- 
ing fluctuations are u EE = £(2)/£(3) = 1.37 for massless bosons as, e.g., 
gluons. Massive bosons have smaller fluctuations with, for example, u n = 1.11 
[50] and u p = 1.01 when T = m n . Massless fermions, e.g. quarks, have 
ujf = 2£(2)/3£(3) ~ 0.91 independent of temperature. 

Resonances are implicitly included in the WNM fluctuations. In the thermal 
limit resonances are found to increase total multiplicity fluctuations [17,24] 
but decrease, e.g., net charged particle fluctuations [51-54]. In high energy 
nuclear collisions, resonance decays such as p — > 2n, uj — > 3n, etc., lead to half 
or more of the pion multiplicity. Only a small fraction r ~ 20 — 30% produce 
two charged particles in a thermal hadron gas [55,51] or in RQMD [56] (see 
also [57]). Not all of the decay particles from the same resonance always fall 
into the NA49 acceptance, 4 < y < 5.5, and the fraction of pairs will be 
smaller; we estimate r ~ 0.1. Including such resonance fluctuations in the BE 
fluctuations gives, similarly to Eq. (10), 

<4*+ R = r 1 -^- + (1 + r)u EE . (18) 
1 + r 

With r ~ 0.1 we obtain uj b t e+r ~ 1.3. In [17] the estimated effect of resonances 
is about twice ours: ujn — 1-5, not including impact parameter fluctuations. 

Fluctuations in the effective collision volume add a further term (N)a(V) / (V) 2 
to uo B I E+R . Assuming that the volume scales with the number of participants, 
oj v /(V) ~ uj Np /(N p ), we find from Eq. (10) that u N = u§ E+R + (n)oj Np ~ 2.1, 
again consistent with the NA49 data. Because of the similarity between the 
magnitudes of the thermal and WNM multiplicity fluctuations, the present 
measurements cannot distinguish between these two limiting pictures. 

3.4 Centrality dependence and degree of thermalization 

It is very unfortunate that the WNM and thermal models predict the same 
multiplicity fluctuations in the NA49 acceptance - and that they agree with 
the experiment. If the numbers from the two models had been different and 
the experimental number in between these two, then one would have had 
quantified the degree of thermalization in relativistic heavy ion collisions. 

The similarity of the fluctuation in the thermal and WNM is, however, a coin- 
cidence at SPS energies. As seen from Fig. (5) the fluctuations in pp collisions 
increase with collision energy and just happen to cross the thermal fluctua- 
tions, ^thermal — 2.2, at SPS energies.^ 

2 As discussed below the thermal fluctuations in positive or negative particles are 
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Fig. 7. The fluctuation in the total number of charged particles (excluding volume 
fluctuations) vs. centrality or energy in the zero degree calorimeter; left are cen- 
tral and right are peripheral nuclear collisions. The curves are linear extrapolations 
between the thermal fluctuations, tothermai — 2.2, in central collisions to pp fluctua- 
tions at SPS, RHIC, and LHC energies, lo pp , expected in peripheral collisions within 
the WNM (see text). 

At RHIC or LHC energies the situation will be much clearer. Here the charged 
particle fluctuations in pp collisions are much larger as seen in Fig. (5), namely 
uj P n = 6.5, 20 at RHIC and LHC energies respectively. The thermal fluctu- 
ations remain as uJthermai — 2.2. Therefore a dramatic reduction in event-by- 
event fluctuations are expected at higher energies at the nuclear collisions 
become more central as shown in Fig. (7). 



This can be exploited to define a "Degree of thermalization" as the measured 
fluctuations at a given centrality relative to those in the thermal and pp limits 

^WNM _ ^exp 

Degree of thermalization = — th ^ rmal , (19) 

uj n — uj n 

which ranges from unity in the thermal limit to zero in the WNM. Whereas 
both u^ NM and u e ^ p may depend on the acceptance the degree of thermal- 

u± — 1.1 in a thermal hadron gas. The fluctuation in total charge is twice that 
due to overall charge neutrality which relates the number of positive to negative 
particles. 
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ization Eq. (19) should not. Contributions from volume or impact parameter 
fluctuations may, however, be centrality dependent and should therefore be 
subtracted. Alternatively, the fluctuations in a ratio, e.g. N_/N + , should be 
taken for limited acceptances. 

At RH1C and LHC it should be straight forward to measure the degree of 
thermalization as function of centrality. This is interesting on its own and 
a necessary requirement for studies of anomalous fluctuations from a phase 
transition. 



3. 5 Enhanced fluctuations in first order phase transitions 

First order phase transitions can lead to rather large fluctuations in physical 
quantities. Thus, detection of enhanced fluctuations, beyond the elementary 
statistical ones considered to this point, could signal the presence of such 
a transition. For example, before it became clear that the chiral symmetry 
restoring phase transition in hot QCD is not a strong first order phase transi- 
tion, it was suggested that matter undergoing a transition from chirally sym- 
metric to broken chiral symmetry could, when expanding, supercool and form 
droplets, resulting in large multiplicity versus rapidity fluctuations [34]. Let 
us imagine that Np droplets fall into the acceptance, each producing n par- 
ticles, i.e., (N) = (N D )(n). The corresponding multiplicity fluctuation is (see 
Appendix B) 

u N = u n + (n)u ND . (20) 

As in Eq. (11), we expect cu n ~ 1. However, unlike the case of participant 
fluctuations, the second term in (20) can lead to huge multiplicity fluctuations 
when only a few droplets fall into the acceptance; in such a case, (n) is large 
and u>n d of order unity. The fluctuations from droplets depends on the total 
number of droplets, the spread in rapidity of particles from a droplet, 8y ~ 



yT/m t , as well as the experimental acceptance in rapidity, Ay. When 5y <C 
Ay and the droplets are binomially distributed in rapidity, ~ 1 — Ay/y tot , 
which can be a significant fraction of unity. 

In the extreme case where none or only one droplet falls into the acceptance 
with equal probability, we have un d = 1/2 and (n) = 2(N). The resulting 
fluctuation is — (N), which is more than two orders of magnitude larger 
than the expected value of order unity as currently measured in NA49. It 
should be said immediately that a much smaller enhancement is realistic as the 
transition probably is at most weakly first order and many effects will smear 
the signal. Yet, this simple example clearly demonstrates the importance of 
event-by-event fluctuations accompanying phase transitions, and illustrates 
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Fig. 8. Qualitative picture of multiplictity fluctuations vs. centrality (total multi- 
plicity or Et)- Anomalous fluctuations appear when a transition to a new state of 
matter (QM) starts at centrality E\ (see text). Curves for different ratios of the 
fluctuations characterizing the two states of matter are shown. 

how monitoring such fluctuations versus centrality becomes a promising sig- 
nal, in the upcoming RHIC experiments, for the onset of a transition. It is the 
hope and expectation that the higher RHIC energies probe deeper into the 
QGP phase by creating higher temperatures and energy densities whereby 
larger regions of QGP are produced. The larger event multiplicities should 
make it possible to improve on statistics and thereby also the ability to detect 
anomalous fluctuations. The potential for large fluctuations (orders of magni- 
tude) from a transition makes it worth looking for at RHIC considering the 
relative simplicity and accuracy (percents) of multiplicity measurements. 

Let us subsequently consider a less extreme model in which a transition leads 
to enhanced fluctuations of some kind. Assume that the total multiplicity 
within the acceptance arises from a normal hadronic background component 
(Nhg) an d from a second component (Nqgp) that has undergone a transition: 

N = N HG + N QGP . (21) 
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Its average is (TV) = (Nhg) + (Nqgp)- Assuming that the multiplicity of each 
of these components is statistically independent, the multiplicity fluctuation 
becomes 



Here, uohg is the standard fluctuation in hadronic matter u>hg — 1- The fluc- 
tuations due to the component that had experienced a phase transition, ujqqp, 
depend on the type and order of the transition, the speed with which the col- 
lision zone goes through the transition, the degree of equilibrium, the subse- 
quent hadronization process, the number of rescatterings between hadroniza- 
tion and freezeout, etc. If thermal and chemical equilibration eliminate all 
signs of the transition, then ujqgp — ojhg- 

The amount of QM and thus (Nq M ) depends on centrality, energy and nuclear 
masses in the collision. For a given centrality the densities vary from zero 
at the periphery of the collision zone to a maximum value at the center. 
Furthermore, the more central the collision the higher energy densities are 
created. The transverse energy, E?, the total multiplicity and/or the energy 
in the zero-degree calorimeter, E Z dc, have been found to be good measures 
of the centrality of the collision at SPS energies. Therefore, it would be very 
interesting to study fluctuations vs. centrality which are proportional to energy 
density. By varying the binning size for centrality one can also remove impact 
parameter fluctuations as discussed above. 

If the energy density in the center of the collision zone exceeds the critical 
energy density for forming QM at a certain centrality, E 1 , then a mixed phase 
of QM and HM is formed. At a higher energy density, where the critical 
energy density plus the latent heat for the transition is exceeded, which we 
shall assume occur at a centrality E 2 , then a pure QM phase is produced 
in the center. These quantities will depend on the amount of stopping at a 
given centrality, the geometry, T c , etc. In the mixed phase E\ < Et < E2, 
the relative amount of QM, (Nqm)/(N), is proportional to both the volume 
of the mixed phase, and the fraction of the volume that is in the QM phase. 
The latter varies in the volume such that it vanishes at HM/QM boundary. 

In Fig. (8) a schematic plot of the fluctuations of Eq. (22) is shown as func- 
tion of centrality for various loqm- Up to centrality E\ the fluctuations are 
unchanged. Above the central overlap zone undergoes the transition to the 
QM/HM mixed phase and fluctuations start to grow when uiqm > ^hm- At 
the higher centrality, E 2 the central overlap zone is in the pure QM phase but 
the maximum fluctuations ooqm are not reached because the surface regions of 
the collision zone is still in the HM phase. On the other hand, if the hadroniza- 
tion of the QM state is smooth and does not lead to enhanced fluctuations 
(i.e. if ujqgp = ujhg), it cannot be observed in such a study. 



LOn = ujhg + (uqgp — uhg) 



(Nqgp) 
(N) 
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The multiplicity fluctuations can be studied for any kind of particles, total or 
ratios. Total multiplicities describe total multiplicities whereas, e.g. the ratio 
7T / (7r + + 7r~) can reveal fluctuations in chiral symmetry. The onset and mag- 
nitude of such fluctuations would reveal the symmetry and other properties 
of the new phase. 



4 Correlations between total and net charge, baryon number or 
strangeness 



By a combined analysis of fluctuations in, e.g., positive, negative, total and 
net charge as well as ratios, the intrinsic and other fluctuations as well as 
correlations can be extracted and exploited to reveal interesting physics as 
will be demonstrated in the following. 



4-1 General analysis of fluctuations and correlations 



Multiplicity fluctuations between various kinds of particles can be strongly 
correlated. As a first example, consider the multiplicities of positive and nega- 
tive pions, N + and AC, in a rapidity interval Ay for any relativistic heavy-ion 
experiment. Similar analyzes can be performed for any two kinds of distin- 
guishable particles. 

The net positive charge from the protons in the colliding nuclei is much smaller 
than the total charge produced in an ultrarelativistic heavy-ion collision. For 
example, (N + ) exceeds (iV_) by only ~ 15% at in Pb+Pb collisions at SPS 
energies. The fluctuations in the number of positive and negative (or neutral) 
pions are also very similar, u N+ ~ oj^_. Charged particle fluctuations have 
been estimated in thermal as well as participant nucleon models [24] including 
effects of resonances, acceptance, and impact parameter fluctuations. By vary- 
ing the acceptance and centrality, the degree of thermalization can actually 
be determined empirically. Detailed analysis indicates that the fluctuations in 
central Pb+Pb collisions at the SPS are thermal whereas peripheral collisions 
are a superposition of pp fluctuations [64]. 

The fluctuations in the total (N ch = N + + N_) and net (Q = N + — N_) charge 
are defined as [54] 

((N + ±N_) 2 ) - (N + ±N_) 2 

(N+ + AC) ~ 

^^ ++ (AW^- ±C ' (23) 
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where the correlation is given by 



°- • (24) 

Fluctuations in positive, negative, total and net charge can be combined to 
yield both the intrinsic fluctuations in the numbers of N± and the correlations 
in their production as well as a consistency check. These quantities can change 
as a consequence of thermalization and a possible phase transition. 

In practice, ~ 0Jn_-> so that the fluctuation in total charge simplifies to 

„ - (N 2 ch ) - (N ch r _ 

UN ch = tt^ u N+ + G , (25) 

and that for the net charge becomes 

U Q = — 7n~\ — = UN + ~ C ■ ^ 



The fluctuation in net charge can be related to the fluctuation in the ratio of 
positive to negative particles 

oo Q ~ (N + /N_)(N ch )u N _ /N+ /4, (27) 

plus volume (or impact parameter) fluctuations [51,53]. The virtue of this ex- 
pression is that volume fluctuations can in principle be extracted empirically. 
Alternatively one can vary the centrality bin size or the acceptance. Further- 
more, the volume fluctuations for net and total charge are proportional to the 
net ((N + — N-)) and total ((N + + iV_)) charge respectively with the same 
prefactor. In the following we shall assume that such "trivial" volume fluctu- 
ations have been removed. 



The analysis has so far been general and Eqs. (23-26) apply to any kind 
of distinguishable particles, e.g. positive and negative particles, pions, kaons, 
baryons, etc. - irrespective of what phase the system may be in, or whether it is 
thermal or not. In the following, we shall consider thermal equilibrium, which 
seems to apply to central collisions between relativistic nuclei, in order to 
reveal possible effects on fluctuations of the presence of a quark-gluon plasma. 



4-2 Charge fluctuations in a thermal hadron gas 



In a thermal hadron gas (HG) as created in relativistic in nuclear collisions, 
pions can be produced either directly or through the decay of heavier reso- 
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nances, p, u;, .... The resulting fluctuation in the measured number of pions 
is 



where f r is the fraction of measured pions produced from the decay of reso- 
nance r, and J2r fr = 1- These mechanisms are assumed to be independent, 
which is valid in a thermal system. 

The heavier resonances such as p°,cu,... decay into pairs of 7r + 7r~ and thus 
lead to a correlation 



Resonances reduce the fluctuations in net charge in a HG in chemical equi- 
librium at temperature T = 170 MeV and baryon chemical potential /if, = 
270 MeV and strangeness chemical potential p s = 74 MeV to ujq = 0.70 
[51,17]. In [52] the value u Q = 0.70 is found. 

In addition, overall charge conservation reduces fluctuations in net charge 
when the acceptance is large and thus increases correlations as will be dis- 
cussed below. 



4-3 Charge fluctuations in a quark- gluon plasma 

A phase transition to the QGP can alter both fluctuations and correlations 
in the production of charged pions. To the extent that these effects are not 
eliminated by subsequent thermalization of the HG, they may remain as ob- 
servable remnants of the QGP phase. As shown in Refs. [52,53], net charge 
fluctuations in a plasma of u, d quarks and gluons are reduced partly due 
to the intrinsically smaller quark charge and partly due to correlations from 
gluons 



where Nf is the number of quark flavors, qf their charges, and N q the number 
of quarks. The total number of charged particles (but not the net charge) can 
be altered by the ultimate hadronization of the QGP. Assuming a pion gas as 
the final state, this effect can be estimated by equating the entropy of all pions 
to the entropy of the quarks and gluons. Since 2/3 of all pions are charged 
and since the entropy per fermion is 7/6 times the entropy per boson in a 



UJ N+ = UJ N _ = f n U„ + f p UJ p + + .... , 




c HG = l -f P + U + 



(29) 




(30) 
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two-flavor QGP 



(N ch )^l((N 9 ) + 7 -(N q )), (31) 

where the number of gluons is (N g ) = (16/9Nf)(N q ) . Inserting this result in 
(30), we see that the resulting fluctuations are ujq = 0.18 in a two-flavor QGP 
(and ujq = 0.12 for three flavors). As pointed out in [53], lattice results give 
ujq ~ 0.25. However, according to [55] a substantial fraction of the pions 
are decay products from the HG, and the entropy of the HG exceed that of a 
pion gas by a factor 1.75 — 1.8. As described in [52] the net charge fluctuations 
should be increased by this factor in the QGP, i.e. ujq ~ 0.33 in a two-flavor 
QGP, whereas it remains similar in the HG, ujq ~ 0.6. 

The above models are all grand canonical ones, i.e. no net charge conservation, 
as opposed to microcanonical models that now will be discussed. If the high 
density phase is initially dominated by gluons with quarks produced only at 
a later stage of the expansion by gluon fusion, the production of positively 
and negatively charged quarks will be strongly correlated on sufficiently small 
rapidity scales. An increased entropy density in the collisions volume will 
lead to enhanced multiplicity as compared to a standard hadronic scenario if 
total entropy is conserved. The associated particle production must conserve 
net charge on large rapidity scales (Ay^l) due to causality because fields 
cannot communicate over large distances and therefore must conserve charge 
within the "event horizon". Therefore the net charge, N c h, is approximately 
conserved whereas the total charge, Q, increase by an amount proportional 
to the additional entropy produced. If the entropy density increases from shg 
to sqgp going from a HG to QGP without additional net charge production, 
fluctuations in net charge will be reduced correspondingly, 

uQGP ^ 8HG_^ (32) 
S QGP 

The resulting fluctuation in net charge is necessarily smaller than that from 
thermal quark production as given by Eq. (30). A similar phenomenon occurs 
in string models where particle production by string breaking and qq pair 
production results in flavor and charge correlations on a small rapidity scale 
[42]. 

If droplets or density fluctuations appear, they are expected not to produce 
additional net charge. Consequently, the net charge fluctuations should still 
vanish ujq ~ whereas oj c h ~2u N+ ~ 2uqgp- 

The strangeness fluctuation in kaons might seem less interesting at first 



3 It is amusing to note that this number gives a very bad (i.e., negative) estimate 
for (N g )/(N q ) in Eq. (31). 
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sight since strangeness is not suppressed in the QGP: The strangeness per 
kaon is unity, and the total number of kaons is equal to the number of strange 
quarks. However, if strange quarks are produced at a late stage in the expan- 
sion of a fluid initially dominated by gluons, the net strangeness will again be 
greatly reduced on sufficiently small rapidity scale. Consequently, fluctuations 
in net /total strangeness would be reduced/enhanced. 

The baryon number fluctuations have been estimated in a thermal model [52] 
in a grand canonical model. It is, however, not known how possible variations 
in baryon stopping event-by-event and subsequent diffusion and annihilation of 
the baryons and antibaryons in the hadronic phase affect these results. If only 
charged particles are detected, but not K°, K°, neutrons and antineutrons, 
the fluctuations have smaller correlations as compared to the total and net 
strangeness or baryon number. 



4-4 Total charge conservation 

Total charge conservation is important when the acceptance Ay is a non- 
negligible fraction of the total rapidity. It reduces the fluctuations in the net 
charge as calculated within the canonical ensemble, Eqs. (28-31). If the total 
positive charge (which is exactly equal to the total negative charge plus the 
incoming nuclear charges) is randomly distributed, the resulting fluctuations 
are smaller than the intrinsic ones by a factor (1 — f acc ), where 



is the acceptance fraction or the probability that a charged particle falls into 
the acceptance Ay assuming full p t coverage. Since charged particle rapidity 
distributions are peaked near midrapidity, charge conservation effectively kills 
fluctuations in the net charge even when Ay is substantially smaller than 
the laboratory rapidity, y^ ~ 6 (11) at SPS (RHIC) energies. Total charge 
conservation also has the effect of increasing uj ch towards 2u)n + according to 
Eq. (25). Similar effects can be seen in photon fluctuations when photons are 
produced in pairs through p° — > 27. In the WA98 experiment, uj 1 ~ 2 is found 
after the elimination of volume fluctuations [49]. 

On the other hand, if the acceptance Ay is too small, particles can diffuse 
in and out of the acceptance during hadronization and freezeout [52]. Fur- 
thermore, correlations due to resonance production will disappear when the 
average separation in rapidity between decay products exceeds the acceptance. 
Each of these effects tends to increase all fluctuations towards Poisson statis- 
tics when Ay^Sy, where 8y denotes the rapidity interval that particles diffuse 




(33) 
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Fig. 9. Acceptance dependence of thermal fluctuations in net charge (uq of Eq. (34), 
lower curves) and total (^N ch , upper curves). Correlations increase from a hadron 
gas (C ~ 0.4), to a QGP (C ~ 0.8) and a gluon plasma (C ~ 1.0) (see text). The 
HG result with a rapidity diffusion of by = 0.8 is also shown for comparison to the 
other curves which use by = 0.5. The large error bar on the NA49 data point is 
not statistical but reflects the uncertainty in the subtraction of impact parameter 
fluctuations from fluctuations in charged particles [11,24]. The corresponding net 
charge fluctuation predicted by UrQMD [58] is shown by open circle. From [54]. 

during hadronization, freezeout and decay. We find approximately 

"^(s^W^sfWK'-*- (34) 

where ujq is the canonical thermal fluctuation of Eqs. (29,30) and ujq V is the 
fluctuation corrected for both 8y and total charge conservation. 

The resulting fluctuations in total and net charge are shown in Fig. (9) as- 
suming oj n+ = oj n ~ 1.1 and 5y = 0.5. As mentioned above, f acc and Ay are 
related by the measured charge particle rapidity distributions [11]. The total 
charge fluctuations in a HG (C = 0.4) from Eq. (31) agree well with NA49 
data [11] after subtraction of residual impact parameter fluctuations. Data on 
charge particle ratios, which do not contain impact parameter fluctuations, 
will be able to test the net charge fluctuations of Eq. (34) to higher accuracy. 
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Predictions from UrQMD are also shown for comparison [58]. The sensitivity 
to diffusion is small as seen in Fig. (9) where for the fluctuations are also 
shown for 8y = 0.8 as recently used in [59]. The curves in Fig. (9) apply to 
RHIC energies as well after scaling 5y with Ay. 



5 Fluctuations in particle ratios 

By taking ratios of particles, e.g. K/tt, it + /it~ , 7i°/ti ± , one conveniently 
removes volume and impact parameter fluctuations to first approximation. 
Simply increasing/decreasing the volume or centrality, the average number of 
particles of both species scales up/down by the same amount and thus cancel 
in the ratio. 



5.1 ratio and entropy production 

Most particles produced in relativistic nuclear collisions are pions and they 
therefore constitute most of the number of positive and negatively charged 
particles. The fluctuations in the ratio and thus the ratio of positive 

and negative particles are intimately related to the fluctuations in net charge 



where uJi p f is the impact parameter or volume fluctuations and ujq are the net 
charge fluctuations as given by Eq. (27). 

The 7r + /7r~ ratio has been studied in detail in [51]. Resonances such as p, u, ... 
decaying into two or three pions correlate the 7r + and 7r~ production as for 
positively and negatively charged particles discussed above. Consequently the 
fluctuation in the ir + /n~ ratio is reduced by ~ 30% in agreement with NA49 
data [11]. 

5.2 K/tt ratio and strangeness enhancement 

To second order in the fluctuations of the numbers of K and tc, we have [24,51] 



[51,53] 
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(N + /N_) u Q + u ipf , 



(35) 



UN-/N+ — 



(N ch ) 



(K/tt) 




(K) 



(Kir) - (K)(tt) 



) 



(36) 
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Fig. 10. Event-by-event fluctuations in the K/ir ratio measured by NA49 in central 
Pb+Pb collisions at the SPS [11]. 

The corresponding fluctuations in (K/n) are given by 

2 _ u K/n u K u n (Kir) - (K)(ir) 

(K/n) (K) + (n) (K)(n) " 1 > 

The fluctuations in the kaon to pion ratio is dominated by the fluctuations in 
the number of kaons alone. The third term in Eq. (37) includes correlations 
between the number of pions and kaons. It contains a negative part from 
volume fluctuations, which removes the volume fluctuations in ujk and u n 
since such fluctuations cancel in any ratio. In the NA49 data [11] shown in 
Fig. (10) the average ratio of charged kaons to charged pions is (K/ir) = 0.18 
and (ir) ~ 200. Excluding volume fluctuations, we take ujk — uj n ~ 1.2 — 1.3 
as discussed above. The first two terms in Eq. (37) then yield D ~ 0.20 — 0.21 
in good agreement with preliminary measurements D = 0.23 [11]. Thus at 
this stage the data gives no evidence for correlated production of K and n, 
as described by the final term in Eq. (37), besides volume fluctuations. The 
similar fluctuations in mixed event analyses D m i xe d = 0.208 [11] confirm this 
conclusion. 

Strangeness enhancement has been observed in relativistic nuclear collisions 
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at the SPS. For example, the number of kaons and therefore also (K/tt) is 
increased by a factor of 2-3 in central Pb+Pb collisions. It would be interesting 
to study the fluctuations in strangeness as well. By varying the acceptance one 
might be able to gauge the degree of thermalization as discussed above. The 
fluctuations in the K/n ratio as function of centrality would in that case 
reveal whether strangeness enhancement is associated with thermalization or 
other mechanisms lie behind. In a plasma of deconfined quarks strangeness is 
increased rapidly by gg — > ss and qq — > ss processes and lead to enhancement 
of total strangeness s + s whereas the net strangeness s — s remains zero. The 
fluctuations in net and total strangeness will qualitatively behave like net and 
total charge, however, with unit strangeness quantum numbers as compared 
to the fractional charges. 

5.3 7T°/7r ± ratio and chiral symmetry restoration. 

Fluctuations in neutral relative to charged pions would be a characteristic 
signal of chiral symmetry restoration in heavy ion collisions. If, during ex- 
pansion and cooling, domains of chiral condensates gets "disoriented" (DCC) 
[36], anomalous fluctuations in tt /^ ratios could results if DCC domains 
are large. For single DCC domain the probability distribution of ratios d = 
7T°/(vr + 7T+ + 7i-) is P(d) = 1/V2d with mean (d) = 1/3 and fluctuation 
ojd = 4/15, i.e. much larger than ordinary fluctuations in such ratios (see Eq. 
(37)) which decrease inversely with the number of pions. 

Neutral pions are much harder to measure than charged pions but with respect 
to fluctuations, it suffices to measure the charged pions only. The anomalous 
fluctuations in ir° due to a DCC are anti-correlated to n ± , i.e. they are of same 
magnitude but opposite sign. A DCC can equally well be searched for in total 
charge fluctuations as in the 7r°/7r ± ratio, except for the troublesome impact 
parameter fluctuations. 

5.4 J/^f multiplicity correlations and absorption mechanisms 

J/ty suppression has been found in relativistic nuclear collisions [60] and it 
is yet unclear how much is due to absorption on participant nucleons and 
produced particles (comovers). Whether "anomalous suppression" is present 
in the data is one of the most discussed signals from a hot and dense phase at 
early times [60] . It was originally suggested that the formation of a quark-gluon 
plasma would destroy the cc bound states [61]. 

In relativistic heavy ion collisions very few J/\l/'s are produced in each colli- 
sion. Of these only 6.9% branch into dimuons that can be detected and so the 
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chance to detect two dimuon pairs in the same event is very small. Therefore, 
it will be correspondingly difficult to measure fluctuations and other higher 
moments of the number of J . 

Another more promising observable is the correlation between the multiplici- 
ties in, e.g., a rapidity interval Ay of a charmonium state ip = J/^fr, ip 1 ', .. (N^) 
and all particles (N) [54]. The correlator (NN^p) — (N)(N^) also enters the 
in the ratio ip/N (see Eq. (36)). The correlator has as good statistics as the 
total number of ip and it may contain some very interesting anti-correlations, 
namely that ip absorption grows with multiplicity N. The physics behind can 
be comover absorption, which grows with comover density, or formation of 
quark-gluon plasma, which may lead to both anomalous suppression of ip and 
large multiplicity in Ay. Contrarily, direct Glauber absorption should not de- 
pend on the multiplicity of produced particles N since it is caused by collisions 
with participating nucleons. 



To quantify this anti-correlation we model the absorption/destruction of ^'s 
by simple Glauber theory 

N± = e -< W = e -jN/(N) ^ (3g) 
tp 

where N§ is the number of J/^'s before comover or anomalous absorption sets 
in but after direct Glauber absorption on participant nucleons. In Glauber the- 
ory the exponent is the absorption cross section times the absorber density and 
path length traversed in matter. The density and therefore also the exponent 
is proportional to the multiplicity iV with coefficient 

= _d\o g N, 

dlogN {N={N) 

In a simple comover absorption model for a system with longitudinal Bjorken 
scaling, it can be calculated approximately [62] 

7-Ef%^W., (40) 



where dN c /dy, a c1 p, v c1 p and tq are the comover rapidity density, absorption 
cross section, relative velocity and formation time respectively. 

On average comover or anomalous absorption is responsible for a suppression 
factor e -7 . It is difficult to determine because only the total ip suppression 
including direct Glauber absorption on participants is measured. 

The anti-correlation is straight forward to calculate when the fluctuations in 
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the exponent are small (i.e. •j^Jun/{N) <C 1). It is 

(NNj,) - (N}(N^) = - 1 uo N {N i) ) . (41) 



It is negative and proportional to the amount of comover and anomalous ab- 
sorption and obviously vanishes when the absorption is independent of mul- 
tiplicity (7 = 0). The anti-correlation can be accurately determined as the 
current accuracy in determining (N^) is a few percent (NA50 minimum bias 
[60]) in each E? bin. 

The anti-correlations in Eq. (41) may seem independent of the rapidity inter- 
val. However, if it is less than the typical relative rapidities between comovers 
and the ifj, the correlations disappear. Preferrably, the rapidity interval should 
be of the order of the typical rapidity fluctuations due to density fluctuations. 

The anticorrelations of Eq. (41) quantify the amount of comover or anomalous 
absorption and can therefore be exploited to distinguish between these and 
direct Glauber absorption mechanisms. In that respect it is similar to the 
elliptic flow parameter for ip [62] for the comover absorption part but differs 
for the anomalous absorption. 



5.5 Photon fluctuations: thermal emission vs. tt° — > 2 7 



WA98 have measured photon and charged particle multiplicities and their fluc- 
tuations versus centrality and binning size. As mentioned above impact 
parameter fluctuations are proportional to the E? binning size; the WA98 
analysis nicely confirms this, and can subsequently remove impact parameter 
fluctuations. The resulting charged particle multiplicity fluctuations with im- 
pact parameter fluctuations subtracted, u>n — (n)ujN p — 11 — 1-2 were shown 
in Fig. (7). 

The fluctuations in photon multiplicities were found to be almost twice as large 
as for charged particles cj 7 — (n)u Np ~ 2.0. This has the simple explanation 
that photons mainly are produced in n° — > 27 decays. The fluctuations are 
then the double of the fluctuations in 7r° to first approximation as seen from 
Eq. (18). 

If the photons were directly produced from a "shining" thermal fireball one 
would expect that they would exhibit Bose-Einstein fluctuations, cj 7 = ui^ E = 
1.37 for massless particles. In addition the 7r°'s in the hadronic background 
will produce photons with uj 1 = uo EE = 2.0. The measured fluctuation in the 
number of photons will therefore lie between these two numbers and can be 
exploited to quantify the amount of thermal photon emission vs. 7T° — > 2 7 
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decay from a hadronic gas 



ivrthermal of) / j ex P 

1 = ZJ_ (42) 

N thermal + N *> 2 .0 - 1.37 ' V ' 

The impact parameter fluctuations must be subtracted from the measured 
photon fluctuations u^ xp by, e.g., taking the ratio of photons to some other 
particle with known behavior. 



6 Transverse momentum fluctuations 



Fluctuations in average transverse momentum were among the first event-by- 
event analyses studied. In a series of papers Mrowczyhski et al. have studied 
transverse momentum fluctuations in heavy-ion collisions with the purpose of 
studying thermalization and other effects. Fluctuations in temperature and 
thus average transverse momentum event-by-event were studied by a number 
of people [15-18] in connection with critical phenoma relevant if the transi- 
tion is close to a critical point. Experimental analyses by NA49 [11,14] reveal 
that a careful evaluation of systematic effects are required before substantial 
equilibration can be claimed in central heavy-ion collisions from transverse 
momentum fluctuations. They also have found strong correlations between 
multiplicity and transverse momentum. 

The total transverse momentum per event 

N 

^ = Eft,*» ( 43 ) 
-i=i 

is very similar to the transverse energy, for which fluctuations have been stud- 
ied extensively [10,8]. The mean transverse momentum and inverse slopes of 
distributions generally increase with centrality or multiplicity. Assuming that 
a = d\og((p t )N)/d\ogN is small, as is the case for pions [63], the average 
transverse momentum per particle for given multiplicity N is to leading order 

(pt) N = (Pt)(l + a(N - (N))/(N)) . (44) 

where (pt) is the average over all events of the single particle transverse mo- 
mentum. With this parametrization, the average total transverse momentum 
per particle in an event obeys (P t /N) = (p f ). When the transverse momen- 
tum is approximately exponentially distributed with inverse slope T in a given 
event, (p t>i ) = 2T, and = 2T 2 = (p t ) 2 /2. 

The total transverse momentum and also the transverse energy contains both 
fluctuations in multiplicity and fluctuations in the individual particle trans- 
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verse momenta and energy (see Appendix C). An interesting quantity is there- 
fore the total transverse momentum per particle, Pt/N, where the multiplicity 
fluctuations are removed to first order although important correlations remain. 

The total transverse momentum per particle in an event has fluctuations 

(N)a(P t /N) = a( m ) + o?(ptfu N + (1 EfoiPtj - ( Pt ) 2 )) . (45) 

The three terms on the right are respectively: 

i) The individual fluctuations c(pt,i) = (p 2 /) ~ (Pt) 2 the main term. In the 
NA49 data, (p t ) = 377 MeV and (N) = 270. From Eq. (45) we thus obtain 
(a(P t /N) 1 ^ 2 1 (p t ) ~ l/y2(N) = 4.3%, which accounts for most of the experi- 
mentally measured fluctuation 4.65% [11]. The data contains no indication of 
intrinsic temperature fluctuations in the collisions. 

ii) Effects of correlations between p t and N, which are suppressed with respect 
to the first term by a factor ~ a 2 . In NA49 the multiplicity of charged particles 
is mainly that of pions for which T ~ (pt)/2 increases little compared with pp 
collisions, and a ~ 0.05 — 0.1. Thus, these correlations are small for the NA49 
data. However, for kaons and protons, a can be an order of magnitude larger 
as their distributions are strongly affected by the flow observed in central 
collisions [63]. 

iii) Correlations between transverse momenta of different particles in the 
same event. In the WNM the momenta of particles originating from the 
same participant are correlated. In Lund string fragmentation, for exam- 
ple, a quark-antiquark pair is produced with the same p t but in opposite 
direction. The average number of pairs of hadrons from the same partic- 
ipant is (n(n — 1)), where n is the number of particles emitted from the 
same participant nucleon, and therefore the latter term in Eq. (45) becomes 
({n(n — 1)) / (n))((p t ,iPt,j^i) ~ (Pt) 2 )- To a good approximation, n is Poisson 
distributed, i.e., (n(n — \))/{n) = (n), equal to 0.77 for the NA49 acceptance, 
so that this latter term becomes ~ ((p tti p t j^i) — (pt) 2 )- The momentum cor- 
relation between two particles from the same participant is expected to be a 
small fraction of cr(p tti ). 

To quantify the effect of rescatterings, the difference between (N) a (Pt/N) and 
(j(pt) has been studied in detail [12] via the quantity 

*(pt)*y/(N)v(P t /N)-yfifaJ. (46) 
As we see from Eq. (45), in the applicable limit that the second and third 
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terms are small, 



1 



(a 2 {pt) 2 uJ N + ((pt,iPtj&) ~ (Pt) 2 )) ■ 



(47) 



In the Fritiof model, based on the WNM with no rescatterings between sec- 
ondaries, one finds <&(pt) — 4.5 MeV. In the thermal limit the correlations in 
Eq. (46) should vanish for classical particles but the interference of identical 
particles (HBT correlations) contributes to these correlations ~ 6.5 MeV [13]; 
they are again slightly reduced by resonances. The NA49 experimental value, 
$(pt) = 5 MeV (corrected for two-track resolution) seems to favor the ther- 
mal limit [11]. Note however that with a ~ 0.05 — 0.1, the second term on the 
right side of Eq. (47) alone leads to $ ~ 1 — 4 MeV, i.e., the same order of 
magnitude. If ((pt,iPt,j&) ~ (Pt) 2 ) is not positive, then one cannot a priori rule 
out that the smallness of $(p<) does not arise from a cancellation of this term 
with a 2 (pt) 2 u>N, rather than from thermalization. 

A comparison of the transverse momentum fluctuations of charged particles 
to those in mixed events, where correlations thus are removed, showed a small 
enhancement of only 0.002 ± 0.002 [11]. It was estimated that Bose effects 
should enhance this ratio by 1-2% but that total energy conservation intro- 
duces an anticorrelation that partially cancels the Bose enhancement [17,18]. 
Experimental problems with two-track resolution have also been estimated 
to lead to a ratio that is 1-2% lower. Consequently, the numbers seem to be 
compatible. 

The covariance matrix between multiplicity and transverse momentum has 
been analyzed by NA49 [11]. Strong but trivial correlations is found due to 
the fact that higher multiplicity gives larger total transverse momentum event- 
by-event. This correlation is removed in the quantity Pt/N and its covariance 
matrix with multiplicity appears diagonal. 



7 Event-by-Event Fluctuations at RHIC 

The theoretical analysis above leads to a qualitative understanding of event- 
by-event fluctuations and speculations on how phase transitions may show up. 
It gives a quantitative description of AGS and SPS data without the need to 
invoke new physics. We shall here look ahead towards RHIC experiments and 
attempt to describe how fluctuations may be searched for. 

General correlators between all particle species should be measured event-by- 
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event, e.g., the ratios [24] 



mm) 



i + 



cu Nj (NjNj) — (Ni)(Nj) 



(48) 



where Nij are the multiplities in acceptances % and j of any particle. Volume 
fluctuations are automatically removed in such ratios, their fluctuations and 
correlations. If the energy deposition, transverse energy or momentum are 
measured, these latter will have additional fluctuation due to the multiplicity 
fluctuations as explained in Appendix C. 

More generally we define the multiplicity correlations between any two bins 



also referred to as the covariance. When i,j refer to two rapidity bins the co- 
variance is also proportional to the rapidity (auto-)correlation function C(yi — 



It is instructive to consider first completely random (uncorrelated or statisti- 
cal) particle emission. For a fixed total multiplicity Nxot, the probability for 
a particle to end up in bin i is pi = (Ni)/N Tot ~ (Ei) / E To t- The distribution 
is a simple multinomial distribution for which 



The % — j result is the well known one for a binomial distribution. The % ^ j 
result is negative because particles in different bins are anti-correlated: more 
(less) particles in one bin leads to less (more) in other bins on average due to 
a fixed total number of particles. 

As shown above there are nonstatistical fluctuations due to various sources: 
Bose-Einstein fluctuations, resonances, etc., and — in particular — density 
fluctuations. As in Eq. (21) we assume that the multiplicity consist of par- 
ticles from a HM and a QM phase. The covariances in Eq. (50) are derived 
analogously to Eq. (22) 

Wij — UJij : HM + {^ij,QM — ^ij,HM> — 7T7T , \OL) 



when (JVj) = (Nj); when different the general formula is a little more compli- 
cated. Now, the hadronic fluctuations oJi^uu is of order unity for i — j, smaller 
for adjacent bins and vanishes or even becomes slightly negative according to 




(49) 



Vj)- 




(50) 
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(50) for bins very different in pseudorapidity or azimuthal angle 0. The QM 
fluctuations can be much larger: iv^qm ~ (^,qm) (see the discussion after Eq. 
(20)). To discriminate the QM fluctuations from the hadronic ones, Eq. (51) 
requires 



The charged particle multiplicity in central Au + Au collisions at RHIC is 
dN ch /drj ~ 500 — 600 per unit pseudo-rapidity [23]. To see a clear increase 
in fluctuations, say Au = ofy- — oj^hm ~ 1, a density fluctuation of only 
{Ni,QM)^j\fN~i — 25 particles are required per unit rapidity corresponding to 
a few percent of the average. By analyzing many events (of the same total 
multiplicity) the accuracy by which fluctuations are measured experimentally 
is greatly improved. Generally, Au ~ 1/ y/N events , and so fluctuations can in 
principle be determined with immense accuracy. 

It may be advantageous to correlate bins with the same pseudorapidity but 
different azimuthal angles since the hadronic correlations between these are 
small whereas QM fluctuations remain. 

No experimental determination of the purely statistical uncertainties associ- 
ated with any one-body distribution — such as multiplicity as a function of 
rapidity — can be performed without measuring and diagonalizing the cor- 
relation matrix dj = (NiNj) — (N,j)(Nj). While it is conventional to assign 
uncertainties according to the diagonal elements M u , the correlations in the 
covariance matrix are required for a correct error analysis and can also reveal 
physical important results. 



8 Summary 

A phase transitions in high energy nuclear collisions, whether it is first order or 
a soft cross-over, density fluctuations may be expected that show up in rapidity 
and multiplicity fluctuations event-by-event. The fluctuations can be enhanced 
significantly in case of droplet formation as compared to that from an ordinary 
hadronic scenario. A combined analyses of, e.g., positive, negative, total and 
net charge, allows one to extract the various fluctuations and correlations 
uniquely. Likewise a number of other observables as charged and neutral pions, 
kaons, photons, J/^f, etc., and their ratios can show anomalous correlations 
and enhancement or suppression of fluctuations. This clearly demonstrates 
the importance of event-by-event fluctuations accompanying phase transitions, 
and illustrates how monitoring such fluctuations versus centrality becomes 
a promising signal, in the upcoming RHIC experiments, for the onset of a 
transition. The potential for enhanced or suppressed fluctuations (orders of 




(52) 
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magnitude) from a transition makes it worth looking for at RHIC considering 
the relative simplicity and accuracy of multiplicity fluctuation measurements. 

An analysis of fluctuations in central Pb+Pb collisions as currently measured 
in NA49 does, however, not show any sign of anomalous fluctuations. Fluc- 
tuations in multiplicity, transverse momentum, K/n and other ratios can be 
explained by standard statistical fluctuation and additional impact parame- 
ter fluctuations, acceptance cuts, resonances, thermal fluctuations, etc. This 
understanding by "standard" physics should be taken as a baseline for future 
studies at RHIC and LHC and searches for anomalous fluctuations and corre- 
lations from phase transitions that may show up in a number of observables. 

By varying the centrality one should be able to determine quantitatively the 
amount of thermalization in relativistic heavy ion collisions as defined in Eq. 
(19) . For peripheral collisions, where only few rescatterings occur, we expect 
the participant model (WNM) to be approximately valid and the degree of 
thermalization to be small. For central collisions, where many rescatterings 
occur among produced particles, we expect to approach the thermal limit and 
the degree of thermalization should be close to 100%. At RHIC and LHC 
energies the fluctuations in the number of charged particles consequently de- 
crease drastically with centrality whereas at SPS energies the two limits are 
accidentally very close. 

Event-by-event physics is an important tool to study thermalization and phase 
transitions through anomalous fluctuations and correlations — as in rain. 
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Appendix A: Damping of initial density fluctuations 



Hydrodynamic flow with Bjorken scaling is stable according to a stability 
analysis carried out in [35] . Bt linearing the hydrodynamic equations in small 
perturbations in entropy density 5s and rapidity 5y around the Bjorken scaling 
solution and looking for solutions in the form of harmonic perturbations, e lkv , 
the hydrodynamic equations could be written in matrix form (Eq. A. 13 in 
[35]) 



5 



5s/ s 





\5y j \-ikc 2 s -(l-c 2 )J\5y 
The eigenvalues of the above matrix 




5s/ s 



(53) 



A ± = ~(l-c2)±^(l-c2)'-c2*', (54) 

always have real negative part for c s fc / and fluctuations are therefore 
damped. For long wave length fluctuations in rapidity and not too soft equa- 
tions of state, c s k > 1 — c 2 , the solution is a damped oscillator. Note that the 
long wave length solution k = reproduces the Bjorken scaling. 

The exact solution for the entropy density fluctuation 

_ = c+e A + ln ( r / r °) + c_e A - hl{T ^ , (55) 
s 

is sensitive to the equation of state through c s , the initial conditions for the 
rapidity density fluctuations (the constants c±), and their wave length k~ l . 

At large times the eigenvalue with the largest real part dominates and 

7 * U ' (56 » 

Here the oscillating factor has been ignored, leaving the power law fall-off of 
fluctuations with exponent 



Min\Re[X ± }\ = 1(1 - c 2 s ) - Re[^\{\ - c 2 ) 2 - c 2 s k 2 } (57) 

One notes that density fluctuations are undamped for soft equation of states 
(c s = 0). They are also undamped if their wave length is long (k ~ 0). 

To estimate the resulting damping we take a typical rapidity fluctuation for 
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a droplet 5y ~ yT/m t ~ 1 discussed above, which corresponds to a wave- 
number k ~ 1. For an ideal equation of state with sound speed c s = 1/V3 
the last term in Eq. (54) is then either imaginary or small and real, and 
the real part of the eigenvalue is dominated by the first term of Eq. (54), 
ite[A±] ~ —1/3. If we take a typical formation time r ~ 1 fm/c and a 
freezeout time t/ ~ 8 fm/c as extracted from HBT studies [1], the resulting 
suppression of a density fluctuation during expansion is a factor ~ = 0.5 
according to Eq. (56). 
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Appendix B: Fluctuations in source models 



As fluctuations for a source model appears again and again (see Eqs. 10,11,18,22) 
we shall derive this simple equation in detail. 

We define the fluctuations for any stochastic variable x as 

u x ^ . (58) 

It is usually of order unity and therefore more convenient than variances. For 
a Poisson distribution, Pjv — e~ a a N /N\, the fluctuation is ojn = 1- For a 
binomial distribution with tossing probability p the fluctuation is un — 1 — P, 
independent of the number of tosses. In heavy ion collisions several processes 
add to fluctuations so that typically uj e ^ p ~ 1 — 2. Correlations can in some 
cases double the fluctuations as, for example, 7r° — > 27 doubles the fluctuations 
in photon multiplicity and net charge conservation doubles the fluctuation in 
total charge. Impact parameter fluctuations further increases the total charge 
fluctuations to uJN ch = 3 — 5 in peripheral nuclear collisions 



Generally, when the multiplicity (N) arise from independent sources (N p ) such 
as participants, resonances, droplets or whatever, 

i=N p 

N= £ rii, (59) 
i=i 

where rii is the number of particles produced in source i. In the absence of 
correlations between N p and n, the average multiplicity is (N) = (N p )(n). 
Here, (..) refer to averaging over each individual (independent) source as well 
as the number of sources. The number of sources vary from event to event 
and average is performed over typically N events ~ 100.000 events as in NA49 
or N events ~ 10 6 in WA98. 

Squaring Eq. (59) assuming that the source emit particles independently, i.e. 
{fiifij) = (rii) (rij) for i 7^ j, the square consists of the diagonal and off-diagonal 
elements: 

(N 2 ) = (N p )(n^ + (N p (N p -l)}(n t }\ (60) 

With (58) we obtain the multiplicity fluctuations 

(N 2 ) - (iV> 2 
UN = (N) = UJn + \ n ) u *t p , 

as in Eq. (10). 
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Appendix C: Fluctuations in the energy deposited 



Many experiments do not measure individual particle tracks or multiplicities 
but instead the energy deposited in arrays of detector segments, E iy in a 
given event. One could also project the energy transversely by weighting with 
the sine of the scattering angle to study fluctuations in transverse energy [7- 
10]. Since particles mostly have relativistic speeds in relativistic heavy-ion 
collisions, the transverse energy is almost the same as the total transverse 
momentum in an event. 



The total energy deposited in the event is 

D 

E Tat = '£E i , ( 61 ) 



and can be used as a measure of the centrality of the collision. The energy 
deposited in each element (or group of elements) is the sum over the number 
of particle tracks (TVj) hitting detector i of the individual ionization energy of 
each particle (e,) 

Ni 

E t = J2tn. (62) 



The average is: (£*) = (Nj)(e). The energy will approximately be gaussian 
distributed, da/dEi oc exp(— (Ei — (Ei)) 2 /2u Ei (Ei)), with fluctuations (see 
Appendix B) 

<** = {E?) { E t \ El)2 =<**+ < e >^ ■ (63) 
Here, the fluctuation in ionization energy per particle 



1, (64) 



depends on the typical particle energies in the detector and the corresponding 
ionization energies for the detector type and thickness. For the BRAHMS 
detectors we estimate u e / (e) ~ 0.3 [65]. This number will, however, depend 
on rapidity since the longitudinal velocity enters the ionization power. As 
these are "trivial" detector parameter, we shall exclude the fluctuations u e in 
most analyses and concentrate on the second term in Eq. (63) which is the 
fluctuations in the number of particles as examined in detail above. 
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